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In this experimental study we show that the coupling between Nb/NbN superconducting mi- 
crowave resonators and their feedline can be made amplitude dependent. We employ this mech- 
anism to tune the resonators into critical coupling condition, a preferable mode of operation for 
a wide range of applications. Moreover we examine the dependence of critical coupling state on 
other parameters such as temperature and magnetic field. Possible novel applications based on such 
nonlinear coupling are briefly discussed. 

PACS numbers: 84.40.Az, 85.25.-j, 42.50.Dv 



Recently, there has been a growing scientific interest 
in nonlinear directional couplers especially in the opti- 
cal regime Ij. It has been shown in a series of publica- 
tions, that operating such devices by means of parametric 
downconversion, may exhibit some important quantum 
phenomena, such as generation of quadrature-squeezed 
states jlj , occurrence of quantum Zeno and anti-Zeno ef- 
fects 0^ non-classical correlations 0, and quantum en- 
tanglement 0. Such effects, if generated and engineered, 
would have immediate and significant implications on a 
wide range of fields, ranging from basic science to opti- 
cal communication systems, and quantum computation 
[^,|^. The dependence of these mentioned effects on op- 
timal coupling strengths and cavity detunings are almost 
evident Q.ll. 

In this study we investigate the coupling between non- 
linear NbN/Nb microwave resonators and their fecdlines, 
as a function of injected power, ambient temperature and 
applied magnetic field. We show that the coupling in 
these resonators could be, to a large extent, dependent on 
externally applied parameters, allowing us thus to tune 
the resonators into critical coupling condition, which is a 
coupling state where the power transferred to the device 
is maximized and no power refiection is present. Such a 
state occurs as both impedances of the device and of the 
feedline match. 

Theoretically critical coupling can be achieved by ei- 
ther tuning the coupling parameters or by altering the in- 
ternal/external losses of the device. Such control was re- 
ported in optical systems , where critical coupling was 
obtained in microring resonator by altering the internal 
loss inside the ring by means of integrated semiconductor 
optical amplifiers. Other critical coupling achieved in a 
fiber ring resonator configuration may be used in opti- 
cal communication systems i.e. switches and modulators 

Our fabricated resonators were assembled in a stan- 
dard stripline geometry, and were housed in a gold plated 
Faraday package made of Oxygen Free High Conductiv- 
ity (OFHC) Cupper. The layout of the Nb and NbN 
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resonators are presented in the insets of Figs01 and |21 re- 
spectively. The resonators were dc-magnetron sputtered, 
at room temperature, on sapphire substrates with dimen- 
sions 34mm X 30mm X 1mm each. A gap of 0.5 mm 
was set between the feedline and the resonators. The 
thickness of the Nb sputtered film is 2200A, whereas the 
thickness of the NbN film is 3000A. Following the de- 
position, the Nb resonator was patterned using standard 
photolithography process and dry etched with photo re- 
sist mask. The sputtering parameters of the NbN su- 
perconducting microwave resonator have been discussed 
elsewhere [lOl |. 

The measurements of the resonators were done using 
the S parameter of a vector network analyzer, while vary- 
ing the input power or the ambient temperature in small 
steps. In Figs. HEl we show a power induced critical 
coupling data obtained at the first mode of the NbN res- 
onator. In Fig. ^we plot the Sn (w) curves of the first 
mode in the complex plane, corresponding to three in- 
put powers 0.944mW, 1.072 mW, 1.14mW in the vicin- 
ity of the resonance, representing p < pc, p — Pc, P > Pc 
states respectively, where Pc is the injected power level at 
which critical coupling is achieved. These three trajecto- 
ries show clearly that the NbN resonator at its first mode 
has transformed from overcoupled (strongly coupled), to 
undercoupled state (loosely coupled) through critically 
coupled condition via input power increase only. In Fig. 
121(a) we show the resonance curves IS"!!! of the three in- 
put powers, presented previously in Fig. ^ in the vicinity 
of the resonance. The resonances were shifted vertically 
for clarity. In Fig. |21 (b) we present the nonlinear evo- 
lution of the I Sill minimum at resonance as the input 
power is increased through the critical coupling power. 
The minimum of the graph ^ —60 dB at Pc = 1.072 mW, 
corresponds to critical coupling condition where there is 
no power refiection. The gradual shift of the resonance 
frequency towards lower frequencies as the input power 
is increased, is plotted in Fig. 121(c). 

In order to identify and characterize the coupling 
mechanism responsible for this dependence on the drive 
amplitude, we apply the following universal expression 
for refiection amplitude of a linear resonator near reso- 
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FIG. 3: Resonator coupling states as defined by their S'ii(aj) 
curves in the vicinity of the resonance loq, (see Eq.Q) cor- 
responding for (a) undercoupled state (b) critically coupled 
state (c) overcoupled state. 



FIG. 1: NbN resonator first mode Sn {uj) curve, drawn in 
the complex plane for three input powers 0.944mW, 1.072 
mW, 1.14mW, with 2MHz span. The resonator changes from 
overcoupled to undercoupled state as the power is increased 
through pc = 1.072 mW. The unit circle and the origin are 
also plotted on the same axis for visual comparison. 
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FIG. 2: Data analysis of NbN resonator first mode, (a) 
|5ii I graphs representing p < pc, p = pc and p > pc cases. 
The resonance curves were shifted vertically by a constant 
offset for clarity, (b) ISn | minimum as a function of input 
power. The minimum of the graph is at pc- (c) Resonance 
frequency vs. input power, (d) Real Sn at resonance vs. 
power difference from critical coupling power. The resonator 
geometry is shown in the inset of subplot (b). 
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5*11(0;) = - 



i{ujQ - w) - (71 - 72) 



i{uja - w) + (71 +72)' 



(1) 



where 71 and 72 are real coupling factors associated 
with the feedline-resonator and the resonator-reservoir 
(dissipation) couplings respectively. However, unlike Un- 
ear oscillators where 71, 72 and ujq are constants, here we 
assume that these factors depend on some externally ap- 
pUed parameter x, such as input power, temperature, and 



magnetic field. This generalization could be justified by 
noting that the \Sii {uj)\ curves measured in the vicinity 
of the resonance, are to a good approximation symmetri- 
cal Lorentzians. The different possible couphng states of 
the resonator as described generally by Eq. are shown 
schematically in Fig. O 

By expanding 71 and 72 near critical coupling drive Xc, 
to lowest order in Sx — x — Xc, we get 71 — j + giSx and 
72 = 7 + ff2<5a^) where 7, gi and 32 are real constants. 

Substituting for 71 and 72 in the vicinity of Xc in Eq. 
at the corresponding resonance uj = ujo, yields 



5*11 (luq) — aiSx + a2 {Sx)^ + O {^5xY 



(2) 



where ai = (52 - 5i) /27, ^2 = (5? ^ sl) /47^- 
In order to evaluate the critical coupling constant 7 
characterizing the resonance presented earlier, we applied 
a fit to the resonance curve at Pc using Eq.(^ and ob- 
tained 7 = 1.3MHz. Moreover, we extracted gi and 52 
by applying a quadratic polynomial fit as in Eq.Q to 
the experimental data of real 5ii at resonance as a func- 
tion of 5p. The data and the quadratic fit are presented 
in Fig. 121 (d), yielding the following nonlinear coupling 
rates gi = -O.lMHz/mW, 52 = 1.2MHz/mW. Having 
52 1 > I51I and g2 > 0, implies that the critical coupling 
in this case, is caused mainly by the increase of the dis- 
sipation coupling as the input power is increased. 

In the Nb resonator, power induced critical coupling 
was attained at 12.589GHz resonance as shown in Fig. 
0] The resonator changes from undercoupled state to 
overcoupled state as the power is increased through 
the critical coupling power at 6.59mW . The coef- 
ficients 7,31,(72 extracted for this case are 48.5MHz, 
1.6MHz/mW, 1.5MHz/mW respectively. Implying sur- 
prisingly that what led to critical coupling in this case, is 
the increase of the feedline-resonator coupling at higher 
rate than the increase of the dissipation coupling rate, 
since we have \gi\ > \g2\, gi > 0. Such increase in the 
feedline-resonator coupling could be attributed to a non- 
linear redistribution of the standing-wave voltage-mode 
inside the resonator relative to the feedline position, lead- 
ing to amplitude dependence of 71 . 
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FIG. 4: Data analysis of Nb resonator critical coupling at 
12.589GHz resonance, (a) |5ii | graphs representing p < pc, 
p = Pc and p > Pc cases. The resonance curves were shifted 
vertically by a constant offset for clarity, (b) 15*11 | minimum 
as a function of input power. The minimum of the graph is 
at Pc ■ (c) Resonance frequency vs. input power, (d) Real 
Sii at resonance vs. power difference 5p from critical cou- 
pling power. The resonator geometry is shown in the inset of 
subplot (b). 



FIG. 5: Data analysis of NbN second mode temperature 
induced critical coupling. (a) 15*11 | graphs representing 
T < TcrityT — Tcrit and T > Tcrit cases. The resonance 
curves were shifted vertically by a constant offset for clarity, 
(b) |Sii I minimum as a function of temperature, the mini- 
mum of the graph is at Tcrit ■ (c) Resonance frequency vs. 
temperature, (d) Real Sn at resonance vs. temperature dif- 
ference ST from critical coupling temperature. 



Based on our experimental demonstration of nonlinear 
coupling between superconducting stripline elements it 
is interesting to examine the feasibility of implementing 
some of the theoretical proposals presented in Rcf. [1-6] 
in the microwave regime. For that end we consider a 
resonator made of two stripline elements having nonlin- 
ear coupling factor g between them. The resonator is 
also coupled to a feed line employed for delivering the 
input and output signals. Consider operating close to a 
resonance having angular frequency ojq and damping rate 
7. A detailed analysis of this model is beyond the scope 
of the present letter. However, employing dimensional- 
ity analysis and some simplifying assumptions one finds 
that the onset of nonlinear instability and bifurcation is 
expected in such a system at input power level of order 
p^^j^ = j/g. In the regime of operation when the in- 
put power is close to p™^ effects such as intermodulation 
gain and quantum squeezing are expected to become 
noticeable IT]. Based on our results with the Nb res- 
onator, we estimate from the experimental values of 7 
and gi (assuming g = gi) that such effects can be imple- 
mented with a moderate input power of order 10~^W. 
On the other hand, effects such as Zeno and anti-Zeno 
require much stronger nonlinear coupling. We esti- 
mate that in the limit of zero temperature such effects 
will become noticeable only when the input power be- 
comes comparable to = ^jfu^j^g^ . Again, assuming 
the experimental values extracted from the data of the 
Nb resonator, one finds p™ = lO^^W - far too high for 
any practical implementation. However, further analysis 
is required to explore other possibilities of implementing 



such effects in the microwave regime. 

In addition to input power tuning, we have succeeded 
to obtain critical coupling using other modes of oper- 
ation, one of which was tuning the ambient tempera- 
ture. To demonstrate this, we applied a constant input 
power to the NbN resonator and measured its second res- 
onance response as we increased the ambient temperature 
in O.OIK steps. The input power apphed was O.lmW 
, lower than of that resonance at 4.2 K, which was 
about 0.209 mW. The resonator changed from over- 
coupled state at T < Tcrit to undercoupled state at 
T > Tcrit- Tcrit was found to be 5.8K. The measure- 
ment results of this mode of operation are summarized 
in Fig. |S1 By applying similar fitting procedures as was 
explained before one can obtain the following constants: 
7 = 4.2MHz, gi = -0.6MHz/K,g2 = 220MHz/K, indi- 
cating that what led to critical coupling condition, in 
this case, was the increase of the dissipation coupling as 
a result of the temperature rise. 

Moreover to better understand the dual relation be- 
tween the input power and the ambient temperature, ad- 
ditional measurement was carried out on the first mode 
of the NbN resonator. We varied the temperature of 
the resonator as a parameter and at each given tempera- 
ture, we scanned the input power in search for the critical 
coupling power pc ■ The measurement result is shown in 
Fig. El where we see that in the case of the first mode 
of the NbN resonator at ~ 2.5GHz, increasing the am- 
bient temperature decreased the input power at which 
critical coupling occurred, and thus divided the T — p*" 
plane into two regions representing the overcoupled and 
undercoupled resonator states, approximately separated 
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FIG. 6: A plot displaying the critical coupling power depen- 
dence of the first mode of the NbN resonator on ambient 
temperature. The critical coupling power decreases as the 
temperature is increased. 



by the guideline shown in the figure. 



Critical coupling condition was also achieved in the 
NbN resonator by means of applied magnetic field. To 
show this we set a constant input power of 1.23niW (the 
critical coupling power at the time of the measurement 
in the absence of magnetic field was 1.259mW), and in- 
creased a perpendicular magnetic field by small steps. 
The results of this mode of operation, unlike power and 
temperature, showed a distinct instability and frequent 
transitions around the critical coupling state as can be 
inferred from Fig. \7\ A possible explanation to this insta- 
ble behavior could be vortex penetration into the grain 
boundaries of our NbN film, as in NbN could be in 
the order of 35mT 
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FIG. 7: Critical coupling achieved by magnetic field at the 
first mode resonance of the NbN resonator while applying 
constant input power of 1.23mW. (a) l^iil vs. frequency 
corresponding to an increased applied magnetic field. The 
resonance curves were shifted vertically by a constant offset 
for clarity, (b) The 15*11 1 value at the resonance frequency. 



In summary, we have succeeded to tune our supercon- 
ducting microwave resonators into critical coupling con- 
dition using input power, ambient temperature and ap- 
plied magnetic field. By data fitting we have been able to 
extract quantitatively the coupling parameters 7,(?i,g2 
and identify the dominant factor responsible for the cou- 
pling change in each case. In addition we have briefly 
discussed some of the possible applications of this vari- 
able coupling mechanism in exhibiting some important 
quantum phenomena in the microwave regime. 
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